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In this work we present a comparison of multiband k • p models, the effective bond-orbital approach, and 
an empirical tight-binding model to calculate the electronic structure for the example of a truncated pyramidal 
GaN/AlN self- assembled quantum dot with a zincblende structure. For the system under consideration, we find 
a very good agreement between the results of the microscopic models and the 8-band k ■ p formalism, in contrast 
to a 6-i-2-band k • p-model, where conduction band and valence band are assumed to be decoupled. This indicates 
a surprisingly strong coupling between conduction and valence band states for the wide band gap materials GaN 
and AIN. Special attention is paid to the possible influence of the weak spin-orbit coupling on the localized 
single-particle wave functions of the investigated structure. 

PACS numbers: 71.15.-m, 73.21.La, 73.22.Dj 



I. INTRODUCTION 

Nitride-based semiconductor nanostructures are promising 
materials due to their potential application in opto-electronic 
and high-power/temperature electronic devices.^ AIN, GaN 
and InN and their ternary and quarternary alloys in principle 
allow the emission of the whole spectrum of visible light from 
red to ultraviolet. Within the past years, increasing research 
interest has been on the investigation of GaN/AlN quantum 
dots (QDs) in order to develop single electron transistors,^ ul- 
traviolet sources^ and detectors.^ 

Group Ill-nitrides can crystallize in the thermodynamically 
stable configuration with a wurtzite crystal structure and in 
the meta- stable modification with a zincblende structure.^ The 
great majority of wurtzite GaN/AlN QDs are grown along the 
polar [0001] direction. These structures exhibit large spon- 
taneous and strain induced polarization. These effects lead 
to a large internal electrostatic field, which is very unique to 
Ill-nitride heterostructures and has a significant effect on the 
electronic and optical properties of QDs. The magnitude of 
the electrostatic built-in field has been estimated to be in the 
order of MV/cm.^-^ Such fields spatially separate the electrons 
and holes, which leads to a reduction of the oscillator strength 
and enhanced radiative lifetimes. ^'^'^ 

In contrast, in the cubic GaN/AlN QD structures, the spon- 
taneous polarization is absent due to the higher crystal sym- 
metry. Furthermore, experimental data indicate that the 
piezoelectric contributions are small. Therefore, GaN/AlN 
QDs with a zincblende structure are expected to have advan- 
tages in optoelectronic devices. Recently, there has been an 
increasing interest in cubic GaN/AlN QDs due to the improve- 
ment of their growth process. ^'^^'^^'^^'^^-^^ In order to under- 
stand the optical properties of cubic GaN/AlN QDs, the in- 
vestigation of the electronic structure of these structures is of 
major importance. For instance, the excitation energies and 
wave functions are crucial ingredients for carrier-carrier^^ and 
carrier-phonon^^ scattering in nitride QD structures. 

Different approaches have been developed to calculate the 
electronic structure of semiconductor QDs. These meth- 



ods range from continuum approaches such as effective 
j^^ggi8,i9,20 . p8,2i,22,23 approximations to atomistic 

models, e.g. tight-binding2i212^ and pseudopotential ap- 
proaches^^'^^'^^. 

The number of available theoretical models for the calcu- 
lation of the electronic structure makes an evaluation of these 
methods with respect to the accuracy of the investigated ma- 
terial properties necessary. 

While in atomistic descriptions the computational effort 
grows with the number of involved atoms, the accuracy of 
continuum models decreases when the structure's characteris- 
tic dimensions reach the length scale of the atomic bonds. On 
the other hand, the continuum models are not limited to a max- 
imum size of the structure. Previously, different k • p-models 
have been compared with an atomistic empirical pseudopoten- 
tial method for different semiconductor sy stems. ^^'^^'^^ These 
investigations revealed various shortcomings in the continuum 
models resulting from the lack of atomistic description or an 
insufficient number of involved bands. However, a compari- 
son of atomistic and continuum models employing the same 
number of involved bands is essential in order to determine 
the accuracy of computationally less demanding models. 

In this work, we perform a careful comparison of vari- 
ous atomistic and continuum methods, namely an empirical 
tight-binding (ETBM) model, an effective bond-orbital model 
(EBOM) and two variations of the k • p approach, to calculate 
the electronic structure of GaN/AlN QDs with a zincblende 
structure. 

Since we use for all methods an equivalent set of input pa- 
rameters, the output of the various approaches can be com- 
pared directly. We focus our attention on the differences in 
the electron and hole wave functions and the corresponding 
single-particle energies. A comparison of optical properties 
such as excitonic absorption or emission spectra, which can be 
obtained in the framework of a full configuration scheme^, is 
beyond the scope of present work. 

In this study special attention is paid to the influence 
of the small spin-orbit coupling on the results, an ef- 
fect which has been commonly neglected in Ill-nitride QD 
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systems. sMsattattittii? However, recent investigations on 
wurtzite Ill-nitride QDs show that neglecting the spin-orbit 
coupHng leads to artificial degeneracies of hole states. ^^-^^ 
The influence of spin-orbit coupling on the properties of 
nanostructures can be expected to be strongly nonlinear in em- 
pirical approaches, as it enters both the common bulk param- 
eter set and the geometry related part of the Hamiltonian. 

The influence of the conduction band - valence band cou- 
pling in the 8 -band k • p-model will be shown to have a sur- 
prisingly large effect on the electron binding energies despite 
the fact that GaN is a wide band gap material. 

This paper is organized as follows: In Sec. [Ill we intro- 
duce the applied methods and their underlying concepts and 
approximations. Section [Till is dedicated to the GaN/AlN QD 
geometry. The following section deals with the electronic 
structure of these systems. The influence of the spin-orbit cou- 
pling will be discussed in detail in Sec. lIVBl In Sec. lIVCl we 
compare the results of the 6- and the 8 -band k • p model. 



II. APPLIED METHODS 

While atomistic models of different sophistication ap- 
proaching the electronic properties of semiconductor nano- 
structures lead to an increasing computational effort with the 
number of involved atoms, continuum models may produce 
strong deviations from results obtained in atomistic simula- 
tions. These deviations are expected to increase with decreas- 
ing characteristic dimensions of the structure. The aim of 
this paper is to provide a comparison of complementary ap- 
proaches to the electronic structure of GaN/AlN QDs in a cu- 
bic structure. 

Previous studies compared the k • p formalism with highly 
accurate but computationally expensive empirical pseudopo- 
tential calculations. ^^'^^ 

In this study we chose two microscopic approaches with 
various levels of approximation which have been constructed 
such that they reproduce the band structure in the for opti- 
cal application relevant region around the F-point. The in- 
vestigated methods are (i) the empirical tight-binding method 
(ETBM), (ii) the effective bond-orbital model (EBOM) and 
(iii) the k • p formalism employing different numbers of bands. 
The choice of equivalent input parameters in the investigated 
methods allows a direct comparison of differences resulting 
purely from the different level of approximation. While the 
ETBM is the most accurate model of the investigated ones, 
the EBOM and the k • p-formalism allow a straight-forward 
study of different material parameters on the electronic struc- 
ture as the input parameter set is fixed for the investigated ma- 
terial system. In this chapter, we will introduce the investi- 
gated methods used to compute the electronic structure of the 
model system. 



A. Empirical tight-binding model (ETBM) 

The key assumption of the tight-binding method is that the 
overlap of atomic orbitals decreases rapidly with the distance 



of their corresponding atoms, i.e., only Hamilton matrix ele- 
ments (TB parameters) between neighboring atoms (typically 
up to the second or third nearest neighbor shell) have to be 
included. For the polar semiconductors GaN and AIN consid- 
ered in this study, the upper valence band is mainly formed 
by the p-orbitals of the anions and the conduction band from 
the s-orbitals of the cations. We therefore apply an Scpl TB 
model,^^ where each anion is described by the outer valence 
orbitals per spin direction: px, Py and p^. The cations are mod- 
eled by a single s orbital per spin direction. Overlap matrix 
elements up to the second nearest neighbors are included in 
our TB model. Following Ref. 41, the spin-orbit component 
of the bulk Hamiltonian H^^^^ couples only ;?-orbitals at the 
same atom. By analytical diagonalization of the TB Hamil- 
tonian f^j. special k directions the electronic dispersion 
is obtained as a function of the TB parameters. Equations 
for the TB parameters can now be deduced in terms of the 
Kohn-Luttinger-parameters (71,72,73), the single particle en- 
ergy gap Eg, the effective electron-mass and the spin-orbit 
splitting Aso at the F-point. Doing so, one TB -parameter has 
to be determined self-consistently to reproduce the L-point en- 
ergy of the split-off band. This parametrization has been veri- 
fied to correctly describe the band structure region around the 
F-point. Since we are dealing with the electronic properties 
of a nanostructure formed from a direct gap semiconductor 
material here, mainly this part of the bulk band structure is of 
importance. Furthermore, due to the large energetic splitting 
(> 1 eV) between the zone center (F-point) and the side val- 
leys (L-, X-point), the quantum confinement will not introduce 
a mixing between these states. Therefore, the region around 
the F-point is expected to dominate the QD states. 

B. Effective bond-orbital model (EBOM) 

In the EBOM the TB orbitals are replaced by so called ef- 
fective orbitals located on the sites of the underlying lattice, 
thus neglecting the atomic basis of the material. With respect 
to the zincblende structure which is considered in this study, 
the underlying symmetry of the original crystal structure is 
changed to that of an fee lattice with effective orbitals 

|I^^^^)eff' i = SiPx^Py^Pz^ (^=T,i, (1) 

on each Bravais lattice site R. We note that this approxima- 
tion gives rise to an artificial change of the symmetry from 
(zincblende) to Oh (fee). 

An advantage of the EBOM approach is that it allows to 
directly relate the underlying TB -parameters with the corre- 
sponding k • p-Hamiltonian. A self-consistent fitting of the 
ETBM parameters to the bulk band structure is therefore not 
needed. 

A first EBOM parametrization by Chang"^^ incorporated 
three-center overlap integrals. This parametrization was re- 
stricted to coupling up to nearest neighbors only, so that 
solely the F-point energies could be fitted to the set of k • p- 
parameters. In the present work, we use the parametrization of 
Loehr^^, which includes hopping up to second-nearest neigh- 
bors to fit the band- structure of the bulk material to the above 
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FIG. 1 : (Color online) Bulk band structure of zincblende GaN in two 
high symmetry directions calculated by the effective-bond-orbital 
method (solid lines), the Scpl tight-binding model (dotted lines) and 
the 8 -band k • p approach (dashed lines). 



mentioned set of parameters. This parametrization addition- 
ally allows for a fit to the X -point energies of the conduction 
band (Xi^), the degenerate hh-/lh-band (X^y) and the split-off 
band (Xsy). The degeneracy of the hh-/lh- band at X is subse- 
quently lifted by the incorporation of spin-orbit-coupling into 
the TB Hamiltonian. 



C. k p-method 

In the k • p-formalism the wave functions are replaced by 
their envelope which no longer resolves individual atoms. By 
this coarse graining treatment it becomes possible to describe 
the electronic structure very efficiently on a continuum scale. 
Within this work, we calculate the electronic structure in a 
basis set of eight complex envelope eigenf unctions. This re- 
sults in an eight-band k • p Hamiltonian^^, given in the Ap- 
pendix. The eight-band formalism can be reduced to two plus 
six bands, if the Kane parameter, £p, describing the coupling 
between the valence band (VB) and the conduction band (CB) 
is set to zero (see Appendix). Taking this parameter into ac- 
count causes modifications of the effective mass and Luttinger 
parameters as well as additional coupling elements U and V 
within the Hamiltonian matrix, H^^^, which are neglected in 
the 6+2-band model. Another common simplification which 
will be checked explicitly in this study is to neglect the spin- 
orbit coupling parameter Aso which is in the order of a few 
meV in the Ill-nitride sy stems. This approximation reduces 
the dimensionality of the Hamiltonian and thus the compu- 
tational effort. The computation of the k • p electron and 
hole wave functions is performed in the plane- wave formal- 
ism within the S/PHI/nX software package. ^^^^ 



TABLE I: Material parameters for zincblende GaN and AlNiSai 



Parameters 


GaN 


AIN 


a [A] 


4.5 


4.38 


£g [eV] 


3.26 


4.9 


A£vb [eV] 


0.8 


0.0 


X{ [eV] 


4.428 


5.346 


[eV] 


-6.294 


-5.388 


[eV] 


-2.459 


-2.315 


-cp Lev J 


25.0 


27.1 


Aso [eV] 


0.017 


0.019 


Me No] 


0.15 


0.25 


71 


2.67 


1.92 


72 


0.75 


0.47 


73 


1.10 


0.85 



D. Comparison of the bulk band structure 

For the parametrization of the bulk band structure, we use 
a set given by Fonoberov et al\ which has been shown to 
accurately reproduce experimental data and recent GqWo- 
calculations^^ around the F-point. 

The parameter sets are given in Tab. Il Figure [T] shows the 
resulting band structure along the L — F — X path for the three 
methods. As can be seen, the agreement around the F-point 
is perfect whereas expected clear deviations arise towards the 
Brillouin-zone boundaries: 

The EBOM reproduces the bulk GaN valence-band structure 
obtained from previous work^^ throughout the Brillouin zone. 
The ETBM results along the F — L direction are in excel- 
lent agreement with the EBOM valence-band structure, while 
slight deviations from the EBOM results are observed at the 
X -point. Remember that the EBOM energies have been fitted 
to this point of the Brillouin zone. 

Furthermore, the EBOM accurately describes the conduc- 
tion band along the T — X direction. In particular, an addi- 
tional maximum along the F — X -direction is reproduced in 
agreement with ab-initio band structure calculations.^^ The 
ETBM conduction band deviates from the EBOM results near 
the X-point since higher conduction-bands are not taken into 
account. 

In contrast to the EBOM and ETBM, respectively, the 8- 
band k • p reproduces the band structure only for small k- 
vectors around F. Similiar to the discussion of the ETBM, 
by taking more bands into account a better agreement of the 
k • p - band structure throughout the Brillouin zone can be 
achieved.— 

However, in the present study we refrain from such an ex- 
tended k • p or TB Hamiltonian to keep the number of in- 
volved bands equal and thus consistent within the investigated 
models. Furthermore, in accordance with the discussion in 
Sec. Ill A[ the F-point character is expected to dominate the 
single particle states in a nanostructure with characteristic di- 
mensions of only a few nm. Indeed, as will be shown later. 
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FIG. 2: Schematic geometry of our model quantum dot. The base 
length b and the height t of the pyramidal frustum determine the dot 
size, while h and d define the size of the supercell. 



even the 8 x 8 k • p Hamiltonian and the Scp\ TB model are 
excellent approximations. 



III. THE MODEL QUANTUM DOT 

Cubic GaN QDs embedded in AIN have been experimen- 
tally investigated by various groups. ^'^^'^^'^^'^^'^^ These stud- 
ies have shown that such QDs grow as truncated pyramidal 
structures. ^^'^^ The QD's are commonly grown in a Stranski- 
Krastanov growth mode, i.e., they form spontaneously when 
the wetting layer (WL) exceeds a critical thickness. 

As pointed out before a major aim of this work is to in- 
vestigate the consistence of the results obtained by the contin- 
uum k • p-approach with the outcomes of (semi-) microscop- 
ical tight-binding approaches explained in Sec. [Ill While the 
accuracy of the continuum k • p model is known to increase 
with the dimensions of the structure, a comparison between 
continuum and atomistic simulations is highly interesting for 
systems with characteristic dimensions in the order of magni- 
tude of the bulk lattice constants. Besides, the small number 
of involved atoms also limits the computational effort of the 
atomistic calculations. 

As a test quantum dot we therefore consider a truncated 
pyramid with a square base length of Z? = 16 a and a height of 
t = Aa only. It is placed on top of a GaN WL with a thick- 
ness of 0.5<2, where a is the AIN lattice constant. This cor- 
responds io t ^ 1.15 nm and b = 1 nm. The dot is oriented 
along the [001] axis. Previous studies of GaN/AlN QDs with 
a wurtzite structure showed no intermixing between Ga and 
Al in these structures.^ ^'^^ Since the structural properties of 
cubic GaN QDs embedded in AIN are found to be similar to 
those of the hexagonal ones^^ we also take in our simulations 
the composition of the nanostructure and the surrounding bar- 
rier material to be pure GaN and AIN, respectively. Only for 
the WL composition experimentally a weak interdiffusion of 
Ga and Al atoms is found.^^ Since the aim of this study is 
on the bound single particle states, which are localized in- 
side the nanostructure, compositional fluctuations in the WL 
region have only minor influence on the bound single parti- 
cle states and will be neglected. Note that the symmetry of 
the outer shape of the QD resulting from the confinement po- 
tential is a C4v-symmetry, while the underlying crystal lattice 



lacks inversion symmetry and reduces the symmetry to C2v 

Within all three approaches, the QD is located in a suffi- 
ciently large supercell to eliminate the influence of the chosen 
supercell boundaries on the single particle states. The conver- 
gence of the eigenstates with respect to the supercell size has 
been carefully checked. 

In the framework of a Scp\ TB model, the C2v symmetry 
of the QD's underlying zincblende structure is naturally in- 
cluded. To set up the Hamiltonian we use the TB parameters 
of the corresponding bulk materials. At the GaN/AlN inter- 
faces we use a linear interpolation of the TB parameters of 
GaN and AIN. Since the nitrogen atoms form a common anion 
lattice the interpolation affects only second nearest neighbor 
elements which are small compared to the nearest neighbor 
contributions. 

The application of EBOM to describe QDs is similar to the 
ETBM approach. The main difference is the restriction to a 
slightly more coarse grained grid where the anion and cation 
positions cannot be resolved. This coarse graining changes 
the underlying lattice from zincblende to fee. Therefore, this 
approach cannot sustain the original C2v symmetry of the 
zincblende structure and increases the number of symmetry 
operations. In this case we are left with a C/^^ symmetry. 

Within the k • p-formalism, the dot and the WL are de- 
scribed by a spatially resolved envelope function. We use the 
bulk GaN parameter set inside the nanostructure and bulk AIN 
parameters for the matrix material. The k • p-formalism does 
not resolve individual atoms and will therefore not reproduce 
the C2v symmetry of the underlying zincblende lattice. The 
QD is simulated on a real-space mesh of 80 x 80 x 80 mesh 
points. 

In this study the focus is on a systematic comparison of 
the different approaches introduced in Sec. HIl rather than on a 
complete description of all aspects of a QD. Therefore, we do 
not consider contributions from strain and electrostatic built- 
in fields in our calculations. Nevertheless, as we will dis- 
cuss in the following section, our results for the single-particle 
level structure are in qualitative agreement with results ob- 
tained in Ref. 8 where strain effects and piezoelectric fields 
have been explicitly taken into account. 



IV. ELECTRONIC PROPERTIES OF CUBIC GaN 
QUANTUM DOTS 

A. Single particle energies and states 

Using the above quantum dot geometry, we applied the 
three approaches to compute the bound hole and electron 
states. In all three approaches, we consistently find a total 
number of eight bound electron states. Due to the large effec- 
tive mass of the hole states and the large valence-band offset 
the number of localized hole states is much higher. While 
higher excited states play an important role in carrier-carrier 
and carrier-phonon interaction, we will restrict our discussion 
on the first four bound electron and hole states since these 
dominate the excitonic and absorption processes in QD struc- 
tures. The energy levels of these states as calculated by the 
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FIG. 3: (Color online) The first four electron and hole single-particle 
energy levels of the GaN QD, as calculated by the effective-bond- 
orbital method (EBOM), the Scpl empirical tight-binding model 
(ETBM) and the 8-band- k • p approach. All energies are given with 
respect to the valence band edge of the AIN. Please note the different 
scaling for the electron and hole energies, respectively. 



three approaches are shown in Fig. [3] 

Fig. [4] shows a top view of the QD geometry and the 
modulus square |v^(r)p of the first four single-particle wave 
functions for electrons and holes, respectively. Each state is 
twofold degenerate due to time reversal symmetry. Qualita- 
tively comparable results have been found by Fonoberov et al. 
for larger truncated pyramidal GaN/AlN zincblende QDsi^ 

Electron states: Comparing the three approaches, the 
single-particle states for the electrons are quantitatively as 
well as qualitatively very similar: Both their energy eigen- 
values and the corresponding symmetry character agree well. 
According to their nodal structure, these states can be classi- 
fied di^ S-, p^-, py- and J-like states, respectively. The electron 
ground state \\fl is ^--like, while the next two states are /7-like. 
In case of the k • p and the EBOM approach the symmetry of 
the system is C4y. Here, the atomic structure of the underlying 
zincblende lattice is not resolved, and the states 1/4 ^^^^ 
energetically degenerate and form linear combinations of the 
form p± = (\/ V^){px + ipy). However, if taking the crystal 
structure into account, as it is done in an empirical TB model 
or a pseudopotential approach, the symmetry is reduced and 
degeneracies are lifted. For the QD considered here, a trun- 
cated pyramidal GaN QD grown along the [001] direction and 
with zincblende structure, the symmetry is C2v This symme- 
try lifts the equivalence between the [110] and [1 10] direction. 
Employing the empirical TB model, we therefore obtain en- 
ergetically non-degenerate px- and py-likc states for i/^l and 
VA3 , respectively. The states i/^ and i/^ are found to be non- 
degenerate with an energy difference of about 0.2 meV. 

This value is much smaller than the energy differences be- 
tween the lower bound states (see Figure [3]). This splitting 



may become more pronounced in other material systems or for 
other QD geometries. Furthermore, inclusion of an atom- 
istic strain field and the piezoelectric potential may also in- 
crease this splitting. 

Hole states: In contrast to the electron states, the hole states 
cannot be easily classified according to their nodal structure. 
This is due to the strong intermixing of the various valence 
bands and prevents a strict classification of the optical selec- 
tion rules on total angular momentum selection rules. This 
finding emphasizes the importance of a multiband approach. 
Qualitatively, we find an excellent agreement of the first four 
hole states, and again the corresponding eigenvalues for all 
three approaches lie within a narrow energy range. However, 
in case of the 8-band k • p approach and the EBOM, the first 
two hole states Xj/^ and i/A^ reveal no spatial anisotropy along 
the [110] and [110] direction, respectively. In the ETBM these 
states show a strong anisotropy along the [110] and [iTO] di- 
rection. This behavior again reflects that only the ETBM ap- 
proach correctly reproduces the C2v symmetry of the system. 
Note that the first two hole states are not degenerate, and all 
three approaches yield a splitting of about 6 meV. This effect 
results (mainly) from the spin-orbit splitting energy Aso and 
will be discussed in more detail in Sec. lIVBl 

The splitting of the first two hole states and the anisotropy 
of these states also has a strong influence on the optical prop- 
erties of these systems. For example, the splitting of the states 
Xf/^ and 1^^^ additional lines in the optical spectra. 

Furthermore, the absence of spatial anisotropy of the states Xf/^ 
and 1/^2 may also lead to a vanishing polarization anisotropy, 
A, for dipole transitions along directions [110] and [110], re- 
spectively. Energy differences in the absolute eigenvalues of 
electron and hole states occur due to the different representa- 
tion of the structure within the investigated models, e.g, within 
the k • p formalism the mesh discretization cannot resolve the 
microscopic representation (i.e. the atomistic nature of the 
interface) of the other two methods. Shifting the QD bound- 
aries slightly by ±3.5 A in the k • p formalism modifies the 
absolute energies of electron and hole states about ±15 meV 
but causes no significant deviations in the energy difference 
of the states with respect to the corresponding ground state. 
Since the electron and hole level structure obtained here from 
the different approaches is found to be similar to the ones in 
Ref . 9 for a wurtzite InN/GaN QD, similar excitonic and emis- 
sion spectra are expected for the present system. Selection 
rules for optical transitions can be analyzed in an analogous 
manner from symmetry aspects, as for example discussed in 
Ref.|9,. 



B. Influence of Spin-Orbit-Splitting 

The spin-orbit coupling has been commonly ne- 
glected in previous studies of nitride based nanostruc- 



turesr 
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since in group Ill-nitrides this contribu- 



tion is of the order of a few meV.^^ The spin-orbit splitting 
is a relativistic effect which increases with the atomic 
number of the atoms. For example for CdSe and Znse this 
splitting is in the order of 0.4 eV.^^ Previous ETBM studies 
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TABLE II: Single-particle energies for the truncated pyramidal GaN 
QD with (Aso ^ 0) and without (Aso = 0) spin-orbit coupling. Each of 
the given states is two-fold degenerate due to spin and time reversal 
symmetry, respectively. 





Aso = 




3-hl-bandk 


P 


4-band k • p 


EBOM 


ETBM 


ei [eV] 


4.5259 




4.4477 


4.3428 


4.4246 


ei [eV] 


4.6768 




4.5759 


4.4621 


4.5677 


^3 [eV] 


4.6768 




4.5759 


4.4621 


4.5677 


eA [eV] 


4.8069 




4.6897 


4.5670 


4.6944 


hi [eV] 


0.6679 




0.6726 


0.7264 


0.7070 


h2 [eV] 


0.6679 




0.6726 


0.7264 


0.7063 


/Z3 [eV] 


0.6622 




0.6627 


0.7193 


0.6962 


/Z4 [eV] 


0.6558 




0.6591 


0.7172 


0.6906 




Aso = 17 meV 




6-i-2-band k 


P 


8 -band k • p 


EBOM 


ETBM 


ei [eV] 


4.5259 




4.4479 


4.3429 


4.3866 


ei [eV] 


4.6768 




4.5761 


4.4623 


4.5152 


^3 [eV] 


4.6768 




4.5761 


4.4623 


4.5154 


^4 [eV] 


4.8069 




4.6900 


4.5672 


4.6284 


hi [eV] 


0.6677 




0.6708 


0.7244 


0.6979 


h2 [eV] 


0.6614 




0.6641 


0.7189 


0.6917 


/Z3 [eV] 


0.6570 




0.6578 


0.7145 


0.6857 


/Z4 [eV] 


0.6505 




0.6539 


0.7123 


0.6812 



for CdSe/ZnSe QDs comparable in shape and size reveal 
a splitting in the first two hole states of several meV.^^'^^ 
Without spin-orbit coupling these states are degenerate in the 
framework of a continuum approach, similar to the results 
reported in Ref. 8 for a truncated pyramidal GaN QD with a 
zincblende structure. Recently, the influence of the spin-orbit 
coupling was shown to break the degeneracy of the first two 
hole states in wurtzite InN/GaN QDs.i22i2^ In this section we 
will therefore discuss the influence of the spin-orbit coupling 
on the electronic properties. 

In Tab. HIl the first four electron and hole energies for 
all three models are given in the presence as well as in the 
absence of spin-orbit coupling. Interestingly, the essential 
Kramers degeneracy resulting in twofold degenerate states left 
aside, all calculations in the present paper bear a lift of degen- 
eracy between the first two hole states of about 5-6 meV for 
the GaN dots, thus in the same order of magnitude as in the 
CdSe dots when spin-orbit coupling is included. 

Since the EBOM parameters are not updated in a self- 
consistent manner, these models provide the opportunity to 
study the influence of the SO-coupling Aso on the single- 
particle states and energies by gradually increasing this pa- 
rameter from zero (which gives the limit of a four-band 
model) up to the final value of 17 meV in GaN and 19 meV in 
AIN, respectively. 

The calculated dependence of the level splitting Is^h^-hi as 
function of Aso is given in Fig. [S] It is linear over the entire 
range. As the bulk band structures of GaN and AIN are hardly 
altered by the relatively small spin-orbit coupling, the main 




Electrons 



EBOM 



8-band k • p 



Holes 




FIG. 4: (Color online) Top view of the truncated pyramidal GaN QD 
structure with the first four bound states for electrons (upper part) 
and holes (lower part). Depicted are isosurfaces of the probability 
density with 10% (red) and 50% (blue) of the maximum value. 
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FIG. 5: (Color online) Energy difference ^h2-hi between the first two 
hole states in dependence of artificially varied spin-orbit coupling as 
calculated with the EBOM and the k • p model. The parameter Ago 
goes over several orders of magnitude, note the logarithmic scale, 
from to an upper value given by Ago = 17 meV for GaN. The level 
splitting depends almost linearly on the bulk spin-orbit-splitting. 



influence on the energy levels stems from the site-diagonal 
incorporation of the spin-orbit coupling into the nanostructure 
Hamiltonian. Additional k • p-calculations give comparable 
results which are also depicted in Fig. [5] 

Hence, we have demonstrated that for the system under 
consideration, despite its relatively small influence on the bulk 
band structure, the inclusion of SO-coupling allows to lift the 
artificial degeneracy of hole states and therefore is of essential 
importance for an accurate description of the single particle 
states. 



C. Influence of CB-VB coupling 

While the eight-band formalism yields good agreement 
with the (semi-) microscopic EBOM and ETBM methods, a 
six-band approach combined with an effective mass model 
provides reliable information for the hole states only. Previ- 
ous studies^ ^'^^ found a strong influence of conduction band- 
valence band-coupling for small and medium band gap mate- 
rials. Despite the fact that GaN has a large band gap of 3.26 
eV, an unexpectedly large coupling is also observed for the 
system considered here: coupling effects between the con- 
duction and valence bands significantly modify the electron 
binding energies. The corresponding values are given in Tab. 
ini However, these couplings have essentially no effect on the 
nodal character of the wave functions. Comparing the differ- 
ence of the 2^^ and 3^^ electron state binding energy to the 
ground state energy, we find 0.1282 eV applying the eight- 
band model, which is in excellent agreement with results from 
ETBM (0.1286 and 0.1288 eV) and EBOM (0.1194 eV) cal- 
culations. The 6+2-band k • p approach which is for the elec- 
trons essentially an effective mass approach, gives an energy 
difference of 0.1509 eV, i.e., it overestimates these energies 
by about 23 meV (17%). The origin of the rather large de- 



viation is due to the value of the Kane matrix parameter, 
which describes the CB-VB mixing effects and contains the 
respective dipole matrix element: Its value is large enough 
{E^ is 25 eV in GaN and 27.1 eV in AIN) to produce non- 
negligible coupling effects even for these wide band gap ma- 
terials. This emphasizes the fact, that not only the energy gap 
is important for a possible decoupling of the conduction- and 
valence-bands, but also the magnitude of the Kane parameter 
is of crucial importance. 

Neglecting the spin-orbit coupling in the 6+2-band ap- 
proach leads to a 3+1 -band model. The energies are given 
in Tab. |II1 While the electron binding energies remain un- 
changed, we find a similar behavior for the hole states as we 
find using the 4-band model by neglecting Aso in the 8 -band 
model: the first two hole states are found to be degenerate 
in this model. When including the spin-orbit coupling, the 
resulting energy difference between the V"^ and the 2^^ hole 
state is 6.3 meV in the 6+2-band and 6.7 meV within the full 
8 -band approach. For the hole states, only small differences 
between the eight- and the six-band model are found. Again, 
we find that the nodal character for the hole state wave func- 
tions as found by the eight-band model is preserved when ne- 
glecting E^. 



V. SUMMARY AND DISCUSSION 

In conclusion, we applied atomistic and continuum mod- 
els, namely k • p-models of different levels of sophistication, 
an EBOM and an Sqp\ ETBM model to derive the electronic 
properties of a zincblende GaN/AlN quantum dot. Starting 
from a set of equivalent parameters fitted to the bulk band 
structure around the F-point for all methods, and applying 
them to the same QD model structure, we find a satisfactory 
agreement between the investigated models for the electron 
and hole wave functions and binding energies. This demon- 
strates that for the GaN/AlN system also the semi-microscopic 
EBOM and the continuum 8 -band k • p model are appropiate 
to describe the electronic properties of nanostructures down 
to feature sizes of a few nm. Small discrepancies between the 
ETBM on the one side and the EBOM and 8-band k • p mod- 
els on the other side are found to be a result of the underlying 
crystal symmetry described correctly in the ETBM only. 
Despite the large band gap of GaN and AIN, we find strong de- 
viations of the electron binding energies between the 8-band 
model and a decoupled 6+2-band approach. These occur due 
to the strong influence of the Kane parameter even in wide 
band gap materials. 

The commonly neglected spin-orbit splitting parameter Aso 
lifts the degeneracy of the first two hole states. Even though 
Aso is much smaller in GaN and AIN than in e.g. CdSe, the 
resulting splitting is in the same order of magnitude in both 
material systems. Neglecting this parameter is therefore not 
suitable in the studied material system. By artificially varying 
Aso, we find a strongly linear correlation between the spin- 
orbit splitting and the energy difference between the first two 
hole states. 
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APPENDIX A: 8-BAND k p HAMILTONIAN 

In a basis set of eight complex wave functions 

^ = (/^, , VA^^^ , , v/[^ , ) 

where r6 denotes the conduction-. Eg the light and heavy hole 
valence band states and Fj denotes the spin-orbit coupling, the 
eight-band k • p Hamiltonian^ can be written as: 
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where the effective mass and the six-band model can be found 
in the 2x2 He for the electron and the 6 x 6 Hy for the hole 
states. Hs denotes the superposition of electron and hole states 
within the eight-band model. The matrix elements are given 
as: 



y 



2mo 



7c {d^ 



R = 



U = 



2mo 



■2^1) 



- /<9v 



^/3 

—i 



Pod, 



P0{d,-idy) 



(A2) 



linger parameters and can be derived from the original Lut- 
tinger parameters if by: 



7c 

72 
73 



= 7^ 



mo_Ep / 2 ^ 1 

Me \Eg ^g + As, 



3£g + As, 
1 



23Eg^l 
1 Bp 



(A3) 



2 3£g + Aso 



Ech and £'vb denote the unstrained conduction and valence 
band offset. Eg = E^h — ^vb is the band gap. Pq is the cou- 
pling parameter between conduction and valence bands, Aso 
denotes the spin-orbit coupling and 



:2m0-| 



(A4) 



Note the minus sign appearing in the S instead of the R ele- 
ment in contrast to Ref . H The Ji denote the modified Lul- 



ls the Kane parameter. 
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